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Abstract 

We continue our investigation concerning the question of whether 
atomic bound states begin to stabilize in the ultra-intense field limit. 
The pulses considered are essentially arbitrary, but we distinguish be- 
tween three situations. First the total classical momentum transfer is 
non- vanishing, second not both the total classical momentum transfer 
and the total classical displacement are vanishing together with the 
requirement that the potential has a finite number of bound states 
and third both the total classical momentum transfer and the total 
classical displacement are vanishing. For the first two cases we rig- 
orously prove, that the ionization probability tends to one when the 
amplitude of the pulse tends to infinity and the pulse shape remains 
fixed. In the third case the limit is strictly smaller than one. This 
case is also related to the high frequency limit considered by Gavrila 
et al. 
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Ionization probabilities of atomic systems in the presence of intense laser 
fields are in general poorly predicted. Intense means here that the field 
intensities are of comparable size in magnitude with the ionization energy of 
the potential and hence conventional perturbation theory ceases to be valid. 
Numerous different methods for theoretical investigations have been carried 
out in order to treat the new intensity regime, such as perturbative methods 
around the Gordon- Volkov solution of the Schrodinger equation [0, ^ ^, 
H, ^ |], fully numerical solutions of the Schrodinger equation [|, [1^, |Tl], [12], |1^, 



|4|, |15|, |T^, Floquet solutions [T^, |1§, |19|, high frequency approximations ||20 



or analogies to classical dynamical systems . Some of these investigations 
have led to the prediction of so-called atomic stabilization, which means 
that the ionization probability is supposed to decrease once a certain critical 
intensity has been surpassed. However, several authors have raised doubts 



28, 


^9, 


3U|. 




29, 


30 



For reviews on the subject we refer the reader to |^, ^ 

In this paper we want to continue our previous investigations 
and answer the question concerning the ionization probability in the limit 
when the field amplitude tends to infinity, while the pulse shape remains 
fixed. Of course strictly speaking one would have to include relativistic effects 
into the analysis at some high intensities and then a proper quantum field 
theoretical treatment is needed. However, the Schrodinger theory with the 
a.c. Stark Hamiltonian is consistent in itself, also in that regime, and in this 
light the limit becomes meaningful. Clearly our analysis does not capture 
the effect of window stabilization, which is the purported phenomenon that 
stabilization only occurs in a certain regime of high intensities and then the 
ionization probability tends to one once this regime is surpassed. 

We consider the Schrodinger equation involving some potential V {x), for 
instance the atomic potential, coupled to a classical linearly polarized electric 
field in the dipole approximation E (t) 

t^^t^ =(^-^ + Vix) + z-E (t)^ i, (f , t) = H {t) i, (f , t) . (1) 

We use atomic units h = e = rrie = ca = 1 and we will mainly adopt the 
notations in [EH]. We now want to state precisely which type of potentials 



and electric fields are included in our analysis. 

Assumptions on V: V (x) is a real measurable function on M'^. To each 
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e > one may decompose V as 



V = Vi + V2 



(2) 



where Vi is in (R^) (i.e. square integrable) with compact support 
and V2 is in L°° (M^) with 



Furthermore we assume that H 
positive bound states. 



= ess sup |V^2(^)| < ^ • 

Ho + V with Ho 



(3) 

4- has no 



Relation (^) means that up to a set of measure zero V2{x) is bounded 
in absolute value by e. We note that the potentials of atoms or molecules 
arising from Coulomb pair interactions belong to this wide class. To obtain 
for instance the decomposition (0) for the Coulomb potential l/\x\ we set 



1 

\x\ 



\x\ 



\x\ 



where Xr{^) is the characteristic function of the ball {x : \x\ < R} of radius 



1 for |x| < -R 
for Ixl > -R . 



Potentials satisfying the above assumptions are Kato small, i.e. for each 
a with < a < 1, there exists a constant (3 = (3{a) > such that 



\\ViP\\<a 



2 ^ 



/3 



(4) 



for all ip e L"^ (M^) with Aip e (M^). The Hamiltonian H is self-adjoint 
on the Hilbertspace (M^) and the domains V (H) and V (Hq) of definition 
of H and Ho agree ||3l|. H is bounded from below and has no positive 
eigenvalues if V decays suitably at infinity 

As for the conditions on the electric field, we assume that it takes on the 
form 

Eit)=Eofit), (5) 
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where / (t) is assumed to be measurable in t with f (t) = unless <t < r. 
We call T the pulse duration, f(t) the pulse shape and Eq the amplitude of 
the pulse E (t) . Further we introduce the quantities 



t t 



bit) = JdsE (s) = EoJdsf (s) = Eobo (t) (6) 



t t 
c (t) = J dsb (s) = tb (t) - j dsE{s)s = EqCq (t) . (7) 



With Cz being the unit vector in 2;-direction, b (r) is the total classical 
momentum transfer and c (r) Cz the total classical displacement. We are now 
in the position to formulate more precisely our assumptions on the electric 
field, that is on the pulse shape f{t). 

Assumptions on E: / (t) is a real measurable non-vanishing function in 
t, with support in the interval [0, r] such that 

bo (r) ^ 0. (8) 

In case the potential possesses a finite number of bound states we only 
assume that 

bo (r)' + Co (rf ^ 0. (9) 

Finally, cq (t) is supposed to be piecewise continuous possibly with a 
finite number of zeros in [0,r]. 

Of course the restrictions of a finite number of bound states excludes the 
Coulomb Potential. However, we would like to remark that in general most 
numerical calculations in this context implicitly also assume a finite number 
of bound states. When projecting on bound states numerically, one is always 
forced to introduce a cut-off. Hence our analysis allows also in that case a 
direct comparison with such computations. The gain in the latter case is 
that when the requirement (§) is relaxed to (P) it allows to include more 
types of pulses such as Gaussian etc. All pulses used in the literature satisfy 
assumption E. 

The ionization probability for any given normalized bound state ip of the 
Hamiltonian H is given by 

V iij) = \\{1-P)U (r, 0) = 1 - \\PU (r, 0) i^f . (10) 
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Here U {t', t) denotes the unitary time evolution operator from time t to time 
t' associated to H{t). Its existence [| follows from results in |3^ (for 
details see 



Further P denotes the orthogonal projection in (M'^) on the space 
spanned by the bound states of H = Ho + V . For more details on the 
precise definition of the ionization probability and its properties we refer the 



reader to ||28|. In what follows f{t) and V will be fixed. 



We now formulate the Main Theorem of this article: 

Theorem 1 With the above assumptions on the electric field E and the po- 
tential V, the ionization probability V for any bound state ip of H = Hq + V 
tends to one for the field amplitude Eq going to infinity 

lim V(tp) = 1 . 

This improves a previous result in |]2S| (see relation (3.31) therein), which 
stated that \im\Eo\^oo'P H^) > 1 — t^c, where c is a constant depending on 
the potential V and on ip only. The proof of this main theorem in case of 
condition (|^) shows that the finite dimensional projector P may in fact be 
chosen arbitrarily. In particular in the case of the Coulomb potential, P may 
be the projector on the space spanned by any finite set of bound states. 
Proof of the Main Theorem: 



To start the proof of the main theorem following [28| we may first rewrite 
the ionization probability as 

V{ip) = l- \\Pexp-ib{T)z-expic{T)p,-U' {T,0)ipf . (11) 

Here U' {f, t) is the unitary time evolution operator associated to the Stark 
Hamiltonian (^ in the Kramers-Henneberger gauge p^ , p5 



H'{t) = -^ + V{x-c{t)e,) . (12) 



^In order to show the existence one actuaUy has to make some additional sufficient 
assumptions on V, namely one assumes that Vi{x — ucz), i = 1,2 are and L°° valued 
continuous functions in u, such that in addition 

W^iu) = ^V,{x-ue,), z = 1,2 
ou 

exists and satisfies ||Wi(u)||p < oo for some 6/5 < p < 4/3 and |jW2(w)||oo < oo uniformly 
in u on compact sets in R. So strictly speaking we have to extend our assumptions on V. 
However, for standard potentials like Coulomb etc. this additional assumption is always 
satisfied and we therefore omitted it above for the sake of clarity. 
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Crucial for the proof of the main theorem will be the next result, that 
in the limit Eq oo the time evolution U' for H'{t) is just the free time 
evolution. We will need this result in the following form 



Theorem 2 For all if e L'^ (R^) 

lim ||(f/'(r,0) -exp-Mo)</'|| =0, (13) 

l-BoHoo 

i.e. U'{t,0) converges strongly to exp(— iri^o) |i^^o| oo. 



The proof of this theorem will proceed in several steps. Before we begin 
with the proof we note that this is essentially Kato's theorem on the strong 



convergence of propagators for time dependent Hamiltonians [^. However 
we cannot use this theorem directly since it is not valid for Hamiltonians with 
Coulomb interaction. Since ||(f/'(r, 0) — exp —itHq) < 2 \\^\\ it suffices to 
prove (|13D for all ip G V{Hq) = V{H), which is a dense set in (R^). 
First we use Du Hamel's formula to write 

{U'{T,0)—exp—iTHo)(p = —i / U' {t, s)V{x — c{s)ez) ey^p —isHo (pds (lA) 

Jo 

with (p G V[Hq). We note that by the spectral theorem exp — isi^o leaves 
T>{Ho) invariant. Therefore from (0) it follows that 

\\{U'{t,0) — exp —irHo) ifW < / \\V{x — c{s)ez) exp —isHo ■ (p\\ ds, (15) 

Jo 

uniformly in Eq. 

To proceed further we use the following 



Lemma 3 For any ip G T> (H) = V (Hq) and all s G [0,r] with co(s) 7^ 
one has 

lim \\V {x- c{s)-fez)(p\\ =0. (16) 

Eo—>oo 
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Proof: It suffices to show that for any ip eT) (H) = V (Hq) 



hm 

|7|^oo 



0. 



We show that for arbitrary small e > the estimate ||V^(x — 76^)^911 < e 
holds for all sufficiently large 7 > 0. 

Since V is Kato small and since —A commutes with translations, the 
potential in the Kramers-Henneberger gauge satisfies a similar estimate 

||^(^-7e.)<^|| < -yV^ +f3M (17) 

with fixed P < 00 and for all 7. 
Indeed, 

II V(f -76^)^11 = \\V{x)exp{i-fp,) ifiW 
< \\Hoexp{i-fp^) ip\\ + P\\exp{ijp^) Lp\\ = \\HoLp\\ + 

where the last equality follows from the fact that Hq commutes with 
the translations. In comparison with we have taken a = 1 and 
chosen P = max(/?(a = 1),1). Hence it suffices to prove ([T6|) on a 
core for Hq which is also a core for H. We recall that C is a core 
for a self-adjoint operator A with domain D {A), if C is contained and 
dense in D (A) with respect to the topology in V (A) given by the norm 
II¥'IId(A) ~ W^fW + llv'll- Indeed, for a given ip G V {Hq) let y?' G C be 
such that 

ll^o(¥^-V^')ll + ll(v^-V^')ll<^. 



Also let 7(£:, if') be such that 



\V {x--fe,)^'\\ < - 



for all 7 > 7(5, if'). 



Then 



IV (f - 762) ip\ 



< \\V (f - 76^) {if - ip')\\ + ||\^ (f - 762) 

< ||iJo(^-^')ll+/5|l(^-^')ll + |<^- 



6 



Now (^^), the set of smooth functions on with compact support, 
is such a core and we will now prove (^) on this core. Assuming that 
G (M^) is normalized, we obtain with the assumptions on V 

||V(x-7e^)v2|| < \\Vi{x-'je^)ip\\ + \\V2{x-'je^)ip\\ 
< \\Vi{x-'je^)Lp\\+e. 

For I7I sufficiently large Vi (x — •yez) ip = and the lemma follows. 



We proceed with the proof of the theorem. Since for fixed ip &V {Hq) the 
map exp(— zsi^o) '■ [0, r] — > "D (Hq) given by s exp(— zsifo)v' is continuous, 
the set S = {exp(— ^si^o)v^}o<s<T compact in V{Hq). Therefore 



\\Vix-cis)e, 







as \Eq\ — > 00 uniformly in ip & S for all s G [0, r] except the finite set where 
co(s) = (see e.g. Now the r.h.s. of (|T5|) for any ip G V^Hq) can be 

bounded by 

||([/'(r,0) -exp-Mo)'/'!! < /" sn^\\V{x- c{s)e,)i)\\ds. (18) 
From (p!7| ) and the definition of S it follows also that 



\\V{x-c{s)e, 
for all G iS and all s G [0, r] with 



A 



< a 



/5||V^I 



for every ip G T>{Hq) uniformly in Eq. 

By the Lebesgue dominated convergence theorem we therefore have that 
the right hand side of (|T8|) tends to zero as \Eq\ — > 00. 

This completes the proof of the theorem. 



Remark 5 From the preceding discussion, it is obvious how to weaken the 
last condition in the assumptions on E. Assume we may divide the in- 
terval [0, r] into 2N + 1 parts as = tq < Ti < • ■ ■ < T2n < 'T2N+1 = t, 
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such that Co(t) vanishes identically in the intervals [T2j, T2j+i] (0 < j < A^) 
and is non-zero except for a finite set in the intervals [t2j+i, T2J+2] 
(0 < J < — 1). Then U' (r, 0) converges strongly to 

(19) 

Since f{t) is by assumption not identically zero, cq {t) is not identically 
zero on [0, r], so U" 7^ exp —irH. 

To prove the main theorem it suffices now by Theorem 2 and the obvious 
estimate 

< \\Pe-'''^^^' ■ e^"(")P^ ■ e-*"^«^|| + \\Pe-'^^^^' ■ e'^^^^"^ ■ (f/'(r, 0) - e"^"^") 

< \\Pe-'''^^^' ■ e'"^^^^^ ■ e-^"^'V|| + || {U\t, 0) 

to show that 



)^l 



lim \\P exp —ib (t) z ■ expici^r) pz ■ exp —irHoipW = . (20) 
|£;o|-*oo 

Here exp —itHq has to be replaced by U" in case remark 5 applies. Since 
exp —itHq leaves V (Hq) invariant, it is enough to show 



lim \\P exp —ib (t) z ■ expic{T) pz(f\\ = . 

l-BoHoo 



(21) 



for aW^eV (H) = V (Hq) in order to prove (ID. 

We now modify some arguments already used in |^ and |2^. First we 
consider the case when bo (r) 7^ 0. Also by assumption we have PH < 0. 

Hence P {H - ^b {rf)' is a well defined operator with norm smaller or 
equal to 2/6 (r)^ . Therefore we have 



I P exp —ib (t) z ■ exp ic (r) p^f 

-1 



P[H-hirf 



1, 



H — -b (r) ) exp ~ib (r) z ■ exp ic (r) pz^p 



< 



biry 



1 2 

H b (r) exp —ib (r) z ■ exp ic (r) pzf 



(22) 
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Inserting the relation 



exp —ic (r) ■ exp ib^r) z ■ H ■ exp —ib (r) z ■ exp ic (r) p^ 

Ho - b (r) p, + h) {rf + V{x-c{t) e,) (23) 



into (|22|) yields 



II P exp —ib (r) z ■ expic (r) Pz^W 
6(r) 

Furthermore 
such that 

Ib.y^ll < II^To^ll + llv^ll . (25) 

Finally we have that || V (x — c (r) e^) <y9|| is uniformly bounded in Eq by 
Lemma 2 and therefore we may control the limit I-EqI ~^ oo in (p^ , i.e. 

the r.h.s. goes as O (^j^^- This concludes the proof of the Main Theorem 

for the case 6o (t) 7^ . 

We now turn to the case when P is a finite dimensional projection and 
^0 (''")^ + Co (t)^ =: ctg 7^ 0- Actually by what has already been proved, it 
would suffice to consider the case boi^r) = 0, Co(r) 7^ only. However, we will 
prove the claim for an arbitrary finite dimensional P not necessarily being 
the projection onto the space spanned by the bound states of H. 

We start with two preliminary considerations. First, by the Campbell- 
Hausdorff formula 

exp —ib (r) z ■ exp ic (r) p^ = exp -b (r) c (r) ■ exp i (c (r) Pz — b (r) z) . (26) 

Now there is always an s such that 

c (r) Pz — b (r) z = EqUq {z cos s + Pz sin s) = EQaoZ{s) . (27) 

Introducing the unitary operator W{s) = exp ^ (p^ + z'^) we may perform a 
Bogoliubov transformation on z 

W{s)zW{s)-^ = Z{s) . (28) 
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Secondly, let v^„(l < n < N) be an ortlionormal basis for the range of P. 
Then 



||Pexp —ih (r) z ■ expic (r) pz(p\ 

N 



n=l 
N 



J] |(l^(s)-Vn,exp^Eoao^ ■ W{s)-^^)\ 

n=X 
N 

E 



n=l 



dx{W{s)~^(p^ (x) {W{s) V) i^) expiEoaoz 



Since {W{s)~^LpJ (x) (ly(s)-V) (^) e (1^^) the right hand side of the last 
equation vanishes in the limit \Eq\ —>■ oo hj the Riemann-Lebesgue theorem, 
which concludes the proof of the Main Theorem. ■ 

Now we turn to the case when 60 (t) = cq (r) = 0. Notice, that if we 
consider linearly polarized light, then for the most common pulse shapes like 
for instance a static envelope, trapezoidal envelope, sine-squared envelope 
etc. the extreme high frequency limit, i.e. u —>■ 00 leads to Bq (r) = 



Co (r) = 0. This limit is needed in order to apply the analysis of Gavrila and 



coworkers |2^, which provides so far the most profound "explanation" for 
the occurrence of stabilization. 

We prove the Second Main Theorem: 

Theorem 4 Let bo (r) = cq (r) = 0. Denote 

p{r)= lim V {^) = - P) exp -irHo ■ ipf . (29) 

\Eo\^oo 

If H has only one bound state (i.e. if P is one dimensional), then p{t) > 
for all T > 0. 

Furthermore p{t) < 1 (at least) for all r G [0, r*], where 



n[{HoiJ,Hoij)-{ilJ,Hoijf] 



21-1/2 
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Proof: We consider the survival probability g(r) = [(■?/', e"*'^^"-?/') |^ and 
observe that p(r) = 1 — g(r) if there is only one bound state. We prove that 
g(r) < 1 for all r > 0. Note that by Schwarz inequahty the bound g(r) < 1 
is trivial. Now 

/"OO 

Jo 

where /i^ is the (nonnegative, absolute continuous) spectral measure associ- 
ated with Hq, 



l\p\<V2\ 

0, A < 0. 



Obviously, J^d^^{X) = 1. It is well known (see e.g. [Q) that |/U^(t)| < 1 
for all r > when the measure /i^ is absolutely continuous. 



The second part of the theorem follows from the estimate of Pfeifer |H3 



We note that due to the Paley- Wiener theorem /I^(t) cannot have com- 
pact support. Therefore the inequality p(r) < 1 must be valid for some 
suitable arbitrary large r > 0. On the other hand, it is well known that 
IKpi^)] ^ Ct~^ for arbitrary N and all large r > 0, since the spectrum of 
Hq is purely transient absolute continuous (see for instance This means 

that in case H has only one bound state, say, the ionization probability p(r) 
will tend to one faster than any power inverse power of r for r ^ oo . 



Example 1 The easiest pulse shape for which theorem 4 applies is f{t) = 
cos(co't), since then co(r = — ) = Bq^t = — ) =0. As a concrete example 
for the potential we choose the Coulomb potential V{x) = —l/\x\. The nor- 
malized wave function of the ground state in the momentum representation 



is given by (see for instance ^5(1) 



P 



,2^2 



such that the survival probability in this case reads 

o 2 



q{r) 



- V 
2 



2\A 



64 



TT 



^2' 



3 

'2' 



IT. 

~2' 
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Here U{a, b; z) denotes a confluent hypergeometric function (see for instance 
[p6| ). So for typical sub-picosecond pulses we obtain for instance g(r = 
400a.M.) = 2.45 10"^ and g(r = lOOOa.u.) = 1.63 10"^. Essential is here to 
note that the survival probability is always non- vanishing and monotonically 
decreasing in r (see figure 1). 

Example 2 We now take the potential to be the attracting point interac- 
tion, often also called the delta potential in three dimensions, (see e.g. Wh ) 
with coupling constant a > 0. This potential has the virtue that it possess 
only one eigenstate 



a e 




-oi\x\ 



2tt \x\ 



with energy 
given by 



[a 



In the momentum representation the wave function is 



a 



TT 



(a^ + p^) 



such that the survival probability turns out to be 



16a2 



dpp'^ 



2 



[a^ + pi- 



,2^2 



1 

vr 



3 1 ira^ 



2'2' 



As figure 2 illustrates, the survival probability decreases monotonically with 
increasing a for fixed pulse duration r. The figure also shows, that for 
increasing r the survival probability decreases. 

We may assume that the Hydrogen atom behaves with respect to the 
energy variation qualitatively the same way as the point interaction. Then 
this example indicates that one should expect that for sufficiently high Ry- 
dberg states the survival probability g(r) will be sufficiently close to 1 even 
for times r 1 ps. 

Conclusion We have investigated the ionization probability in the ex- 
treme intensity limit for three different situations. The first analysis presumes 
that the classical momentum transfer &o(''") is non- vanishing and allows essen- 
tially all common potentials. Since the condition 6o(''") 7^ excludes a wide 
range of possible pulses, we also studied separately a situation for which we 
only demand that not both the classical momentum transfer hoij) and the 
total classical momentum transfer cq{t) vanish simultaneously. In addition 
we have to demand for this case that potential only possess a finite number of 
bound states. This is similar to the situation in many numerical calculations. 
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in which one is also forced to introduce a cut-off at some level when projecting 
onto bound states. In both cases we find that the ionization probability V for 
any bound state of the Hamiltonian H = Hq + V for the field amplitude Eq 
going to infinity is going to one. This excludes in our opinion the possibility 
of stabilization for these situations, apart from window stabilization. 

Finally, we considered the situation in which hQ{T) = Cq{t) = and find 
indeed the possibility of stabilization. For the most common pulses, which 
involve linearly polarized light, this case corresponds to the high frequency 
limit of Gavrila and coworkers [^. We conclude that our analysis is consis- 
tent with the "high frequency picture" and that stabilization is only to be 
expected in this latter case. 
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Figure 1: Survival probability g(r(a.M.)) after the time r for the ground state 
of the Hydrogen atom under the free time evolution. 



Figure 2: Survival probability qa{T{a.u.)) for fixed pulse duration r = 200a. m. 
dotted hue, r = AOOa.u. dashed hne, r = 1000a. sohd hue, for the bound 
state of the three dimensional delta potential under the free time evolution 
as a function of the coupling a. 
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